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Abstract This paper describes a particularly transparent derivation of the 
Hawking effect for massive particles in black holes. The calculations are per- 
formed with the help of Painleve-Gullstrand's coordinates which are associated 
with a radially free-falling observer that starts at rest from infinity. It is shown 
that if the energy per unit rest mass, e, is assumed to be related to the the 
Killing constant, k, by k 2 = 2e — 1, then e, must be greater than i. For par- 
ticles that are confined below the event horizon (EH), k is negative. In the 
quantum creation of particle pairs at the EH with k = 1, the time component 
of the particle's four velocity that lies below the EH is compatible only with the 
time component of an outgoing particle above the EH, i.e, the outside particle 
cannot fall back on the black hole. Energy conservation requires that the par- 
ticles inside, and outside the EH have the same value of e, and be created at 
equal distances from the EH, (1 — ri n = r out — 1). Global energy conservations 
forces then the mass of the particle below the EH to be negative, and equal to 
minus the mass the particle above the EH, i.e, the black hole looses energy as a 
consequence of pair production. 

1. The Metric 

A transparent derivation of Painleve-Gullstrand's coordinates (used in this 
paper) is the following: Consider a radially, free-falling observer that starts at 
rest from infinity. For this observer, the equations for r as a function of the 
proper time r, and the equation for the Schwarzschild time, t, as a function of 
r, are respectively (cf. Hartle Eqs. (9.38), (9.40)), 




where r* and t* are two constants. The equations are written in geometrized 
units. It follows from Eq. (l.a) that the second term in Eq.(l.b) is equal to 
t — r*. This relation suggests the following transformation of coordinates, 



the constants in Eqs. (l.a) and (l.b) have been ignored because only dt is 
important. From Eq. (2) we obtain, 

dt = dr-dr^—j (1--) . (3.a) 
which can also be derived directly from the equations, 
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valid for a radial plunge with no kinetic energy at infinity. The new metric can 
be written, 



ds' 



= -(l- dr 2 + 2 (^pj 2 drdr + dr 2 + r 2 (d9 2 + sin 2 6dcf> 2 ) (4) 



The radially ingoing and outgoing light rays are found to be (dr/dr; it follows 
from the geodesic equations that r is an affine parameter), 
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The condition ds 2 — —dr 2 requires that, 

{^)^ + 2^ffdrdr + dr 2 ^, (6) 

the value of dr/dr (= —w) obtained from Eq. (6) agrees with Eq. (l.a). 
2. Particle Orbits 

Expression for a particle's energy per unit rest mass. 

We introduce the following change of notation: the time coordinate, r, in 
Eq. (4) for the metric will be designated hereafter by t, whereas r will be a 
particle's proper time. For radial motions the equations for u r = dr/dr and 
u l = dt/dr are, 



u • kv = - l) u< + WU r = -k, (u r f = k 2 + ^ - 1, W = (^) ^ (7) 

where u = (u',?/*), kv is the Killing vector, namely (1, 0), and A: is a constant 

(the Killing constant). For metrics that are independent of t, an integral of the 
geodesic equation exists, that is given by the scalar product of (1, 0) with u, 
which is Eq. (7.1), namely the first equation in (7). Equation (7.2) follows then 
from the normalization condition for u, i.e., u 2 = — 1. In the derivation of this 
equation, terms in the product ku r appear, but they cancel. In Eq. (7.1), w is 
dimensionless and it is clear that one can assign to k the dimension of a velocity 
(if dimensions are introduced the first term would be multiplied by c.) We 
assume that k and e, the energy per unit rest mass, are related by, k 2 = 2e — 1, 
and find that Eq. (7.2), i.e., the second equation in (7), can be written, 
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e = 1 + i y--T- (8) 

Little need is there to praise Eq. (8)! Because k is constant it follows from 
k 2 = 2e — 1, that e is also a constant, and if k = 1, also e = 1. 

Equation (8) is an expression for the conservation of the particle's energy as 
it moves along a time-independent space-time geometry. 

Particle orbits for some important values of k. 

In this section we study particle orbits only for values of k = 0,1, because 
then (unlike other cases as, e.g., for 1 > k > 0, in particular) solutions can 
easily be found. If k = 0, and for radial motions, the equation for u is found to 
be, 

u= {to (w 2 - , - (w 2 - 1)* ,0,o| ,fc = (9) 

We expect that for physically meaningful solutions the particle's proper time 

increases with the coordinate time, and have therefore discarded the solution 

_ i_ 

with dt/dr = —w (w 2 — l) 2 . Equation (9) is valid only below the event horizon 
and it follows from Eqs. (7.2) and (8) that the value of — u r given by Eq. (9) 
is the minimum fall velocity, and furthermore that e = |. Therefore inside 
the event horizon, e must be larger than ^. Below the event horizon, the 
potential energy decreases, but this decrease is compensated by the increase 
in the minimum allowed value for u r . 

We proceed now to calculate the orbits for |fc| = 1. The value k = 1, 
corresponds to the orbit of a free-falling observer that starts at rest from infinity 
The value of e for all the solutions with |fc| = 1 is unity. It will be shown that 
k = — 1, defines the orbit of a particle, with the same energy per unit mass as 
the free falling observer, but confined below the event horizon. It is convenient 
to use an orthogonal system of coordinates associated with the falling observer. 
It can be readily verified that the vectors, 

e = {l,-w}, ex = {0,1}, (10) 
form an orthonormal basis; the first one being timelike. 

Notice that in this basis, the vector u in Eq. (9), for k = 0, and the Killing 
vector, kv, can be written, 

u = — — ; , kv = e + wei = (l,w) (11) 
(w 2 - l) 2 

where we have adopted the following convention: {v ,vi}, and (v ,vi), denote 
vectors in the coordinate and orthonormal basis respectively. It is clear that in 
Eq. (11), u 2 = -1, and that u-kv= —k = 0. 

Returning to the case \k\ = 1, the equations that need to be satisfied are, 

u • kv = —k, - (u 1 ) 2 + {u r f = -1, |fc| = 1. (12) 
The four solutions to Eqs. (12) are, 
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u=(l,0), k = l; u=(-l,0), k = -l (13.6) 




We do not expect the particle's proper time to decrease while the coordi- 
nate time increases, we are therefore left with only three physically meaningful 
solutions. In the coordinate basis the components of these three solutions are, 



(14.o) 
(14.6) 



The solution with k = 1, in Eq. (13. b), shows that u is the velocity of a 
particle that is stationary with respect to the falling observer, as expected from 
our choice of basic vectors in Eq. (10). The same conclusion can be reached in 
the coordinate basis (cf. the first equation in (14. a)), because it follows from 
Eq. (l.a), that dr/dr = —w. 

Equation (14. b) represents a sinking particle with e = 1 and has no phys- 
ical meaning above the event horizon because there, dt/dr < 0. Conversely, 
the second equation in (14. a) represents an outgoing particle with no physical 
meaning below the event horizon. Here, the minimum descent velocity is equal 
to (w 2 — l) 2 from Eq. (9), and w in Eq. (14. b) must satisfy the inequality 
w > (w 2 — l) 2 , which is indeed the case. 

3. Particle Pair Creation and Hawking Effect 

The Killing vector, kv = (1,0), is timelike above the event horizon (out), 
and spacelike below (in). In the quantum creation of a pair of particles, energy 
conservation requires that (Hartle, page 291), 

-u in ■ kv - u out ■ kv = k in + k ut = 0. (15) 
Here, u in , u out are the four- velocities of the created pair. Above the event 

horizon, k outl must be positive because it is proportional to the particle's en- 
ergy measured by an observer with velocity kv. Therefore ki n — — u^ n • kv 
must be negative and equal to —k out . For values of \k\ = 1, it follows from 
Eq. (14. b), that the velocity of the created particle below the event horizon 
must be, u in = {(w 2 n + l) / {w 2 n — l) , — Win} ■ From the time dependence of the 
solutions in Eq. (14. a), it is apparent that u out must then be taken equal to 
{— {w 2 ut + l) / (w 2 ut — l) ,w ou t} , which is the outgoing solution for k = 1. The 
particle above the event horizon cannot fall back into the black hole. Because 
Win and w out are both very approximately equal to one, it is straightforward 
to show that u\ n — u* Mt requires that Eq. (16.1) below, be satisfied. Equation 
(16.2) follows from the relation, e = (k 2 + l)/2, 

where Ti n and r out are the radial coordinates of the particles forming the pair, 
an intuitively attractive result. 
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In Eq. (16.2), ej„ and e ou t are the energies of the respective particles divided 
by their mass. Because e;„ = e ou t, global energy conservation requires that the 
mass of the particle below the EH be equal to minus the mass of the particle 
above the EH, i.e., the energy of the particle below the event horizon must be 
negative in agreement with Schutz and Carlip's interpretation of the Hawking 
effect. The particle with negative mass survives for a finite amount of time 
before reaching the center of the black hole. But it is an unobservable particle 
and provides the formalism with the necessary degrees of freedom that allows for 
the correct interpretation of an observed particle at infinity escaping from the 
black hole. It would of course be of great interest to understand what happens 
for values of k such that < k < 1 because then the outside particle cannot 
escape to infinity. 

We calculate now the ratio dt/dr, where t is the Schwarzschild time, and 
t is the proper time of the particles at the event horizon. From Eq. (3b) it 
follows that dt/ dtt me t r i c \=l/ (l — iu 2 ) , and then from Eq. (14a) we obtain for 
an outside particle, 

dt_ = (1 + w 2 out ) = 1 

dr (l-wt ut f 2(l-w out f 

In a theory of particle creation, the proper time should play the relevant 
role. Assume then that at the event horizon, N particle-pairs are produced 
in a time dr. Equation (17) shows that for the outside observer, N particles 
will have been produced in the incomparably larger time, dt — dr/2 (1 — w ou t) , 
which suggests a weak observed productions of particles. However, the Hawking 
radiation from a black hole is also very weak. Field theory calculations show 
that black holes emits as though it were a black body with temperature, 
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kbT =WGMY (18) 

the notation being standard. The temperature, T, is truly the physical tem- 
perature of the black hole, not merely a quantity paying a role mathematically 
analogous to temperature in the laws of black hole mechanics (cf. Wald, p. 12). 
Mini black holes excepted, the emission is weak, and the creation of parti- 
cles with finite mass, even neutrinos, must be weaker still. However, as the 
temperature increases, during the final stages of evaporation, the creation of 
particle-pairs with finite mass, could conceivably, become important. It is clear 
however that the answer to this issue lies far beyond the scope of this paper, 
and can only be obtained with the help of field theories capable of calculating 
particle creation in a curved space time (see, e.g., Parker). 

4. Acknowledgements 

Enlightening comments by Professor James Hartle on the first version of 
this article are acknowledged. I am grateful to Drs. Bertrand Chavineau, Roger 
Clark, and Robert Low for helpful discussions. 

5. References 



Carlip, S.: Re: Hawking radiation and vacuum fluctuation, http://sci.tech-archive.net/sci.physcsresearch/200' 
Hartle, J.B.: Gravity, Addison Wesley. California, (2003) 
Parker, L.: Quantized Fields and Particle Creation in expanding Universes, 
Phys.Rev. 183, 1057 (1969) 



5 



Schutz, B.: A First Course in General Relativity, Cambridge University 
Press (2007) 

Wald, R.M.: The Thermodynamics of Black Holes, Living Rev. Relativity A, 
(2001), 6. URL (cited on 1,24,2008) [http://www.livingreviews.org/lrr-200"T6| 



6 



